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Abstract 

A variety of dimensionless numbers related to momentum and heat transfer are useful in 
CVD analysis. These numbers are not traditionally calculated by directly using reactor 
operating parameters, such as temperature and pressure. In this paper, these numbers have been 
expressed in a form that explicitly shows their dependence upon the carrier gas, reactor 
geometry, and reactor operation conditions. These expressions were derived for both 
monatomic and diatomic gases using estimation techniques for viscosity, thermal conductivity, 
and heat capacity. Values calculated from these expressions compared well to previously 
published values. These expressions provide a relatively quick method for predicting changes 
in the flow patterns resulting from changes in the reactor operating conditions. 
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1. INTRODUCTION 


Chemical vapor deposition (CVD) is an industrially important process involving both 
homogeneous and heterogeneous chemical reactions in which a solid film is deposited from the 
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gas phase onto a solid substrate. Since the often complex flow patterns inside a CVD reactor 
oan greatly affect the purity and uniformity of the deposited fllm, it is important to be able to 
predict how they will change in response to changes in reactor conditions. While a complete 
computational model is required to provide detailed information regarding a process as 
complex as CVD, dimensionless numbers related to die flow cm provide useful information on 
the relative magnitudes of various effects, such as how flow patterns vary with reactor operatmg 
conditions. These numbers are typically expressed in terms of viscosity, thermal conductivity, 
heat capacity, and density, as shown in Table 1. These common formulations do not readily 


show their dependence on common operating parametem, such as temperature and pressure. 

Using estimation methods based on kinetic theory for viscosity, thermal conductivity, and 
heat capacity, expressions have been derived for some common dimensionless numbers useful 
in describing the flow in a CVD reactor, die Prandd number (Pr), the Knudsen number (Kn), 


the Reynolds number (Re), thermal Peclet number (Pe). thermal Grashof number (Gr), mid the 
thermal Rayleigh number (Ra). Each expressions is given in form: 


dimensionless number - AxBxCxD 


where: A = a constant dependent on the units used for other parameters in the 

expression 

B = parameters related to the carrier gas used 

C = parameters related to reactor geometry 

D = parameters related to reactor operation conditions 
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Using these expressions, it should be possible to quickly determine the likely effects on the 
flow patterns of altering the reactor operation conditions. It is important to remember that the 
expressions for thermophysical properties and dimensionless numbers used in this paper are 
valid only when all variables are expressed in the units defined m this paper. 


2. ESTIMATION OF THERMOPHYSICAL PROPERTIES 


Since the chemically reacting species in CVD reactors are often present in very low 
concentrations, a dilute solution is assumed in this analysis so that the thermophysical 
properties of the gas in the reactor are those of the carrier gas. Ideal gas behavior is also 
assumed, which is a valid assumption at high temperatures and low pressures. CVD reactors 
typically operate at or below atmospheric pressures and at elevated deposition surface 
temperatures, thus satisfying the conditions for the ideal gas assumption. 

The Chapman-Enskog theory may be employed to estimate the Newtonian viscosity if 
it can be validly assumed that: 1) the system is comprised of a dilute gas where only binary and 
elastic collisions occur; 2) the collisions are described by classical mechanics; and 3) the 
intermolecular potential is spherically symmetric [2], These assumptions apply to low 
pressure, high temperature monatomic gases, but have also been found to be quite good for 
polyatomic gases [1]. From this theory, Newtonian viscosity may be expressed as:[2] 


(I = (2.669 x 10 


-6. (MT) 05 

) 9 


( 1 ) 


where: = Newtonian viscosity (kg/m-s) 

M = molecular weight (kg/kgmoles) 
T = temperature (K) 
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ct = collision diameter (angstroms) 

£2^ = collision integral (dimensionless) 


k T 

The collision integral is a function of the dimensionless temperature, , as shown in Figure 
1 131: 


«. . iff) 

where: k B = Boltzmann constant (1.38049x10 J/K) 

T= temperature (K) 

e = characteristic energy of interaction between molecules (J) 

Excellent curve-fits to the data in Figure 1 are given by the following formulae [4]: 


In Q,* = 0.557 - 0.145 ln(^) - 0.405|erf( S^l if 

knT 

-5-^6 

8 

(3) 

/k B T\ -01487 

= U615 (-r) 

!^Z>6 

8 

(4) 

Expressions for estimating the thermal conductivity of a gas must take into account 
whether the gas is monatomic or polyatomic. For a monatomic gas, the Chapman-Enskog 
theory may be used to obtain an expression for the thermal conductivity [2]: 

8.322x10 f =-=J 

k — - — (monatomic gas) 

a 2 ^ 

(5) 

where: k= thermal conductivity (W/ m-K) 

M = molecular weight (kg/kgmoles) 
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T = temperature (K) 
a = collision diameter (angstroms) 

Q = collision integral (dimensionless) 
r 

For a polyatomic gas, the thermal conductivity may be estimated by using the modified Eucken 
correlation [2]: 


k - 32 ^^^ - 1 j + 1.77^ (polyatomic gas) (6) 

where: k = thermal conductivity (W/m-K) 

(i = Newtonian viscosity (kg/m-s) 

R = universal gas constant (83 14 J/kgmole-K) 

M = molecular weight (kg/kgmoles) 

C p = constant pressure heat capacity (J/kg-K) 

While the heat capacity tends to be nearly constant for monatomic gases, it may vary 
with temperature for polyatomic gases and can be expressed in terms of a polynomial m 
temperature. Such an expression, however, would result in expressions too complex for the 
purposes of this paper. When ideal gas behavior is assumed, simple expressions may be used 
for the heat capacity of monatomic and diatomic gases [5]. 


5R 

(monatomic gas) 

(7) 

2M 



7R 

(diatomic gas) 

(8) 

2M 




where: C p = constant pressure heat capacity (J/kg-K) 

R = universal gas constant (8314 J/kgmole-K) 
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M = molecular weight (kg/kgmoles) 


Using equations (1), (5), (6), (7), and (8), the viscosity, thermal conductivity, and heat 
capacity for argon, helium, hydrogen, and nitrogen were estimated at the temperatures 280 K, 
500 K, 1000 K, 1500 K, and 2000 K. The estimated values were compared with published 
values [5-7] in Figures 2-5. The estimated values for argon appear to be quite good, being 
within 5%, 7%, and 0.1% for viscosity, thermal conductivity, and heat capacity, respectively. 
For helium, the estimated viscosity is within 6% of the published value. The estimated thermal 
conductivity is within 3% of the published value at 280 K. Data points for the thermal 
conductivity at 1500 K, 2000 K, and 2500 K had to be extrapolated. The difference between 
the curves is about 22% at 2500 K; this difference arises both from the estimation technique 
and the extrapolation, but the extent of the contribution of each one is unknown. The published 
values for the heat capacity are those given by equation (7), so the estimated and published 

values are the same. 

For hydrogen, the difference between the viscosity curves is less than 11%, less than 
24% for the thermal conductivity curves, and less than 12% for the heat capacity curves. 
Extrapolated values were used for viscosity and heat capacity at 2500 K, and for thermal 
conductivity at 2000 K and 2500 K. For nitrogen, the viscosity curves differ by less than 6%, 
the thermal conductivity curves by less than 24%, and the heat capacity curves by less than 
19%. Extrapolated values were used for viscosity and heat capacity at 2500 K. The differences 
between the curves for both hydrogen and nitrogen result from both the estimation and the 
extrapolation. Again, the individual contributions of each is unknown. 

Equation (6) for the thermal conductivity, known as the modified Eucken correlation, 
usually predicts values which are higher than those observed experimentally [2], While the 
estimated values for the four gases examined here are higher than published values at lower 
temperatures, they are lower than the published values at the higher temperatures. 
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Since the purpose of this paper is to be able to get a rough idea of the type of flow inside 
a CVD reactor at various conditions, the errors in the estimated values for the thermophysical 
properties are low enough to allow a useful analysis. 


3. EXPRESSIONS FOR DIMENSIONLESS NUMBERS 


The Prandtl number is the ratio of the Newtonian fluid properties [8]: 


Pr _ v _ (9) 

= a k 

Substituting (1), (5), and (7) into (9) yields a constant for the Prandtl number: 


Pr = 0.667 


(monatomic gas) (10) 


Substituting (1), (6), and (8) into (9) yields a constant for the Prandtl number: 


Pr = 0.690 


(diatomic gas) (11) 


The Knudsen number is used to determine if a continuum approach is valid in analyzing 
a system. Generally, continuum fluid dynamics can be used to study flow for Kn<0.01. For 
0.01<Kn<10, the flow is transitional. For Kn>10, a molecular approach must be used to 
describe the flow [9]. The Knudsen number is expressed as: 



( 12 ) 


where: A = mean free path (m) 
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L = characteristic length (m) 


In a dilute gas mixture in which species B is the predominant gas, any of the other 
species wiU most likely collide with a molecule of species B. The distance a molecule of A 
travels before colliding with molecule B may be expressed as [10]: 

_ 4V / M A + M B y °- 5 (13) 

X ” 7tN B (a A + a B ) 2 ^ M b 

where: V = volume (nr) 

M a = molecular weight of species A 
M b = molecular weight of species B 
N b = number of molecules of species B 
a, = collision diameter (angstroms) 

Using the equation of state for an ideal gas, the quantity ^ ma y be expressed as: 

Nb = x b pa o (14) 

V RT 

where: N B = number of molecules of species B 

V = volume (m^) 

x B = mole fraction of species B 
P = pressure (Pascals) 

Aq = Avogadro’s number (6.023xl0 +26 /kgmole) 

R = universal gas constant (8314 J/kgmole-K) 

T = temperature (K) 

Substituting equations (13), and (14) into equation (12) yields: 
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( 15 ) 


Kn = 


4R rMA^v°- 5 __L_ri¥-i- > i 
nA 0 l M b J ( 0/t + a B ) 2 ^ x B P J 


By multiplying out the constants, with appropriate corrections for using the more 
common units of Angstroms for a, instead of meters, the final expression for the Knudsen 

number is obtained: 


Kn = 1.758x10 


-3 


7 m a +m b) 

- 0.5 1 


fJL ' ) 

V m b J 

* (°A + °B ) 2 - 



(16) 


The Reynolds number is the controlling parameter for forced convection. In pipe flow, 
for instance, the flow is always laminar for a Reynolds number less than 2000 [11]. Forced 
convection can be neglected if (Re Pr) 05 « 1 [1 1]. The Reynolds number is given by: 

Re = U7) 

V 

where: (v> = average velocity (m/s) 

L = characteristic length (m) 
v = kinematic viscosity (m /s) 

The kin ematic viscosity may be expressed as [1]. 

„ = B < 18 > 

p 

where: p = Newtonian viscosity (kg/m-s) 

p = density (kg/m 3 ) 
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From the equation of state for an ideal gas, density may be expressed as: 

_ PM (19) 

p RT 

where: p = density (kg/m ) 

P = pressure (Pa) 

M = molecular weight (kg/kgmole) 

R = universal gas constant (8314 J/ kgmoles-K) 

T = temperature (K) 

The average velocity may be obtained by summing all velocities over a cross-section, 
then div iding by the cross-sectional area [1], Such an expression would be too complex for the 
purposes of this paper. For a reactor with a constant cross-sectional area, the mass flow rate 

into the reactor for steady flow may be expressed as: 

( 20 ) 

m = p(v>A cs 

where: m = mass flow rate (kg/s) 

p = density (kg/m 3 ) 

(v) = average velocity (m/s) 

= cross-sectional area (m ) 

By substituting equation (19) into equation (20) and rearranging the terms, and 
approximation for the average velocity is obtained. 

. . mRT (2D 

<v> = PMA^s 


10 



where: <v) = average velocity (m/s) 

m = mass flow rate (kg/s) 

R = universal gas constant (8314 J/ kgmoles-K) 

T = temperature (K) 

P = pressure (Pa) 

M = molecular weight (kg/kgmole) 

\ cs = cross-sectional area (m ) 

In order to obtain the average velocity, the mass Sow rate, m , should be the average 
mass Sow rate through the reactor. Since many CVD reactors use a dilute solution of reactive 
species in a carrier gas, the mass Sow into the reactor is essentially equal to the mass Sow out 
of the reactor. In this case, the average mass Sow rate should be essentially equal to the inlet 
mass Sow rate, which can be obtained from the settings of the mass Sow controllers for the 

reactor. 


Substituting equations (1), (18),(19), and (21) into equation (17) yields: 


Re = 3.747x10" 


(c\) 

f L\ 

f * \ 

m 

M° 5 j 

[a cs ) 

~0.5 

v 1 / 


( 22 ) 


The Peclet number is the thermal equivalent of the Reynolds number. Forced convective 
heat transfer can be neglected if Pe° 5 * 1 [8]. The Peclet number may be expressed as: 

Pe = (23) 
a 

where: (v) = average velocity (m/s) 
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L = characteristic length (m) 
a = thermal diffusivity (m /s) 


The thermal diffusivity is defined as: 


a = 


k 

pCp 


( 24 ) 


where: k = thermal conductivity (W / m-K) 

p = density (kg/m 3 ) 

C P = constant pressure heat capacity (J/kg-K) 


An expression for the Peolet number for a monstomio gas may be obtained by 
substituting equations (5). (7), and (19) into (24), then this result, along with (21), into (23): 


'o 2 Q^ 

Pe = 2.498xl0'’| — of 
M 


L 

A C 


( • \ 

m 


r.0.5 


(monatomic gas) 


(25) 


Similarly, an expression for the Peolet number for a diatomic gas may be obtained by 
substituting equations (6), (8), and (19) into (24), then this result, along with (21), into (23): 


Pe = 2.586x10“ 


a 2 ^ 


M 


0.5 


L 
A r 


m 

rraO.5 


(diatomic gas) 


(26) 


The Grashof number, which may be interpreted as the ratio of the momentum flux by 
free convection to that by diffusion, is useful in determining the magnitude of free convective 
effects. The Grashof number plays the same role in free convection as the Reynolds number 
does in forced convection [12). No turbulence is expected for a Grashof number well below 
lxlO* 5 [131- The ratio Gr/Re" is often used to describe the occurence of buoyancy effects. 
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ii 


The value of k can depend on the reactor configuration and the value of the Reynolds number. 
The Grashof number may be expressed as [1]: 

gPjL AT (27) 

Gr = 5 — 

v 

where: g = gravitational constant (9.8 m/s ) 

(3 X = coefficient of volume expansion (K l ) 

L = characteristic length (m) 

AT = characteristic temperature difference (K) 
v = kinematic viscosity (m /s) 

(3 t is defined as: 



Assuming an ideal gas, Pj becomes equal to 1/T. 


By substituting equations (1), (18), (19), and (28) into equation (27), the following equation is 
obtained: 


Gr = 2.031xl0 3 g(Ma 4 Q |Jl 2 )(L 3 ) 



(29) 


When g = 9.8 m/s, equation (29) becomes: 


Gr = 1 . 9903 xl 0 4 (Ma 4 Q |l 2 )(L 3 )| 


ATP 

T 


(30) 
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The Rayleigh number, which may be interpreted as the ratio of the thermal flux by free 
convection to that by diffusion, characterizes the convective state of an enclosed gas subjected 
to a destabilizing temperature gradient, such as a heated bottom and cold upper wall [14]. This 
is often the situation in horizontal CVD reactors with walls actively-cooled by water and a 
heated deposition surface. Natural convective heat transfer can be neglected if 
(Ra • Pr)° 25 « 1 [8]. Longitudinal convective rolls have been found to appear in horizontal 
CVD reactors for Ra >1708. The transition from laminar to turbulent flow occurs at Ra greater 
than about 2300 [15]. The Rayleigh number may be expressed as: 

„ gP T L 3 AT ( 31 ) 

Ra — — — — — 
va 

where: g = gravitational constant (m/s ) 

P T = coefficient of volume expansion (K _1 ) 

L = characteristic length (m) 

AT = characteristic temperature difference (K) 

v = kinematic viscosity (m /s) 

* 

a = thermal diffusivity (m /s) 

By substituting equations (1). (5), (7), (18). (19). (24). and (28) into equation (31), an 
expression of the Rayleigh number for a monatomic gas is obtained: 

/ 2 \ 

Ra = 1 . 354 xl 0 3 g(Ma 4 nJ)(L ! ) (monatomic gas) (32) 

V T J 

When g = 9.8 m/s, equation (32) becomes: 

/ 2 \ 

Ra = 1.3267xl0 4 (Mo‘V)(L 3 ) (monatomic gas) (33) 

V T J 
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By substituting equations (1). (6), (8), (18), (19), (24), and (28) into equation (31), an 
expression for the Rayleigh number for a diatomic gas is obtained: 

Ra = 1.402xl0 3 g(Mo 4 nj|)(L 3 )f-— £-1 (diatomic gas) (34) 

l T ; 

When g = 9.8 m/s, equation (34) becomes: 

/ 2 \ 

Ra = 1 . 3739 xl 0 4 (Mo 4 £i ll 2 )(L 3 ) (diatomic gas) (35) 

l T / 


A comparison of equations (17) and (23) shows that the Reynolds number and Peciet 
number differ only in their denominators; comparing equations (27) and (31) shows that the 
Grashof and Rayleigh numbers also differ only in their denominators. Since the momentum 
diffusivity, v , and the energy diffusely, a , are always of the same order of magnitude for 
gases [11], it is expected that substitution of expressions for these quantities into expressions 
for the dimensionless numbers would yield essentially identical expressions for the Reynolds 
number and Peciet number (as a comparison of equation (22) with equations (25) and (26) 
shows), as well as essentially identical expressions for the Grashof and Rayleigh numbers (as 
a comparison of equation (29) with equations (32) and (34) shows). 

All the dimensionless numbers examined in this paper have temperature as a variable 
in their expressions. The temperature varies widely in a CVD reactor, making it difficult to 
determine tire temperature to use in calculating the numbers. Since tire material is deposited 
on a substrate surface, flow patterns over this deposition surface will often be of the most 
interest. As an approximation, tire temperature used to calculate the dimensionless numbers 
may be taken as the average between tire reactor wall and the substrate surface: 

^sub + ^wall (36) 
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where: T sub = substrate surface temperature (K) 

T wa u = reactor wall temperature (K) 


The expressions for the Grashof and Rayleigh numbers contain a temperature 
difference. Again, since the main area of interest in likely to be over the deposition surface, the 
temperature difference can be taken to be the difference in temperature between the reactor wall 

and deposition surface: 

^ (37) 

AT = T sub -T wall 

where: T sub = substrate surface temperature (K) 

I'll = reactor wall temperature (K) 


5. SAMPLE CALCULATIONS USING ESTIMATION EXPRESSIONS 

As an example of the use of the various estimation expressions, consider a very simple 
horizontal CVD reactor with a square cross-section, 0.05 meters on a side. The top and side 
walls are maintained at a temperature of 280 K. while the bottom wall, where the substrate 
surface is located, is maintained at a temperature of 1800 K. The gas in die reactor consists of 
an helium earner gas with a mole fraction of silane equal to 3x10^. The operating pressure is 
1.01325xl0 f5 Pa. The mass flow rate into the reactor is equal to 2.3x10 kg/s. 

To begin, values for various parameters needed to calculate the dimensionless numbers 

will be obtained from the given conditions: 
x b = 1-3x10 4 = 9.997x10* 1 
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Acs = 0.05 x 0.05 = 2.5xl0' 3 m2 
From (36): T = (280 + 1800)/2 = 1040 K 
From (37): AT = 1800-280 = 1520 K 

From (21): <v> = (2.3xl0' 6 )(8314)(1040)/(1.01xl0- 5 )(4.003)(2.5xl0' 3 ) = 2x10 2 m/s 

k _Bj _ ^040 / 10.2 = 102; therefore, use (4) to calculate £2^ : 

8 

= 1.1615(102)"° 1487 = 0.584 

Except for the Knudsen number, which uses physical parameter values for both helium 
and cii.~- the physical parameters of helium are used in all calculations. Since the silane is 
present in such a low concentration, the properties of the gas in the reactor are essentially those 

of helium. 

In calculating the Knudsen number, let species A be silane and species B be helium. 
Physical parameters for helium may be found in Table 3. For silane [16]: 


M = 32.1 kg/kgmole 
ct = 4.084 angstroms 
e/k B = 207.6 K 

Equation (16) is used to calculate the Knudsen number: Kn = 3x10 . This value 

implies that continuum fluid dynamics can be used to study the flow. 

The Reynolds number is calculated from equation (22): Re = 1. This value implies 
that the flow in the reactor wiU be laminar. (Re ■ Pr) 03 = 0.8 , which indicates that forced 

convection cannot be neglected in this problem. 

The Peclet number is calculated from equation (25): Pe = 7. This value indicates that 

heat transfer by convection and diffusion are about the same order of magnitude. This is 
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confirmed by Pe" = 0.8 , which shows that forced convection heat transfer cannot be 
ignored. 

The Grashof number is calculated from equation (30): Gr = 2000. Tbrbulent flow is not 

expected for this value of the Grashof number. 

The Rayleigh number is calculated from equation (33): to - 1300. No longitudinal 

convective rolls should occur: however, this value is getting close enough to the critical value 
of Ra = 1708 that it is not inconceivable that such rolls would be present. The equations used 
to calculate these numbers contain a number of simplifying assumptions, and cannot be 
expected to yield exact results. (R. Pr) 025 = 6. which implies tha, free convective heat 

transfer cannot be ignored. 


4. EVALUATION OF EXPRESSIONS FOR DIMENSIONLESS NUMBERS 

Values for dimensionless numbers calculated using the expressions in this paper were 
verified against previously published values that were calculated using thermophysical 
properties obtained from a computational model of a CVD reactor [17). The comparison is 
shown in Table 2. Molecular weights and values for the collision diameter used to compute 
the dim-"»iQ"'ess numbers are shown in Table 3. The values computed using the expressions 
in this paper agreed well with those computed from values obtained from the computational 
model, such as thermophysical properties, average temperature, and average velocity. Since 
both the computational model and the expressions above calculated these parameters using 
estimation expressions based on kinetic theory, their agreement serves as a check on the 
derivation of the expressions for dimensionless numbeis in terms of process parameters. 


5. CONCLUSIONS 
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Dimensionless numbers useful for analyzing flow in CVD systems have been expressed 
in terms of parameters commonly under the control of the system operator, such as temperature 
and pressure. These expressions clearly indicate the dependence of these numbers on system 
parameters, and allow a quick determination of the flow regime inside the reactor that is 
sufficiently accurate for many applications. These expressions also allow a quick 
deter mination of the sensitivity of the flow to changes in system parameters. 
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notation 


Symbol 

Acs 

Ao 

C P 

8 

k 

L 

1 ^ 

Ni 

P 

R 

T 

AT 

V 

a 

Py 

e 

K 

P 

P 

a 

v 


Description 

cross-sectional area 

Avogadro’s number 

constant pressure heat capacity 

gravitational constant 

thermal conductivity 

Boltzmann constant 

characteristic length 

molecular weight 

number of molecules of species i 

pressure 

universal gas constant 
temperature 

characteristic temperature difference 
volume 

mole fraction of species i 
thermal diffusivity 
coefficient of volume expansion 

characteristic energy of interaction between molecules 
exponent in ratio Gr/Re 
Newtonian viscosity 
density 

collision diameter 
kinematic viscosity 

collision integral for viscosity and thermal conductivity 


Units 

m 

kgmole ' 1 

J/kg-K 

m/s 2 

W/m-K 

J/K 

m 

kg/kgmoles 

molecules 

Pa 

J/kgmoles-K 

K 

K 

m 

dimensionless 

m 2 /s 

K 4 

J 

dimensionless 

kg/m-s 

kg/m 3 

angstroms 

m 2 /s 

dimensionless 
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Table 1: Dimensionless numbers useful in CVD flow analysis 


Name 

Prandtl 

Knudsen 

Reynolds 

Peclet 

(thermal) 

Grashof 

(thermal) 

Rayleigh 

(thermal) 


Description Physical Interpretation 



momentum diffusivity 
thermal diffusivity 



mean free path 
characteristic length 



momentum flux by convection 
momentum flux by diffusion 



thermal flux by convection 
thermal flux by diffusion 


gP T L 3 AT 


momentum flux by free convection 
momentum flux by diffusion 


„ gM-’AT 

Ra = 

va 


thermal flux by free convection 
thermal flux by diffusion 
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Table 2: Comparison of dimensionless numbers 


Dimensionless 

Number 

Equation 

Used 

Estimated 
(P=2.7 kPa) 

Estimated 
(P=101 kPa) 

Model [17] 
(P=2.7 kPa) 

Model [17] 
(P=101 kPa) 

Pr 

(ID 

6.9x1 0 _1 

6.9x1 O' 1 

1 

1 

Kn 

(16) 

2 xlO 4 

6 xlO" 6 

lxlO -4 

lxlO -6 

Re 

(21) 

2 

2 

3 

3 

Pe 

(25) 

2 

2 

1 

1 

Gr 

(28) 

0.2 

200 

0.2 

200 

Ra 

(33) 

0.1 

150 

0.2 

200 
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Table 3: Physical parameters for gases [1] 


Gas 

M (kg/kgmoles) 

a (Angstroms) 

e/k B 

argon 

39.944 

3.418 

124.0 

helium 

4.003 

2.576 

10.2 

hydrogen 

2.016 

2.915 . 

38.0 

nitrogen 

28.02 

3.681 

91.5 
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Figure 1: Collision integral for viscosity and thermal conductivity calculations as afunction of 
dimensionless temperature for the range: a) 0 < -|- < 10 ; b) 10 < — < 400 - 

Figure 2: Comparison of estimated and published property values for argon [5-7]: a) viscosity; 
b) thermal conductivity; c) heat capacity. 

Figure 3: Comparison of estimated and published property values for helium [5-7]: a) viscosity; 
b) thermal conductivity; c) heat capacity. The estimated and published values for 
helium heat capacity are the same. 

Figure 4: Comparison of estimated and published property values for hydrogen [5-7]: 
a) viscosity; b) thermal conductivity; c) heat capacity. 

Figure 5: Comparison of estimated and published property values for nitrogen, 
a) viscosity; b) thermal conductivity; c) heat capacity. 
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